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Abstract.  Our  turbulent  lengthscale,  velocity  and  diffusivity  scalings  are 
compared  with  data  from  other  numerical,  laboratory  and  field  experiments, 
Comparison  is  made  with  reference  to  the  turbulence  intensity  measure 
e/uN2.  We  showed  that  our  turbulent  lengthscale  and  velocity  results  are 
consistent  with  measurements  from  the  experiments  considered,  and  that 
the  rms  turbulent  lengthscale  Lt  is  independent  of  the  rate  of  dissipation 
of  turbulent  kinetic  energy  when  e/vN2  > 300.  A diffusivity  modelled  in 
terms  of  an  advective  buoyancy  flux,  6,  is  found  to  reproduce  our  direct 
measurements  of  Kp  in  the  experiments  considered  when  e/uN2  > 300  . It 
is  shown  that  modelling  Kp  as  0 .2e/N2  is  a poor  parameterisation  of  the 
advective  buoyancy  flux  model  in  all  the  experiments  considered,  and  that 
at  large  e/uN2,  this  parameterisation  can  over-predict  the  true  Kp  by  two 
orders  of  magnitude.  This  overprediction  is  discussed  in  terms  of  a mixing 
efficiency  and  it  is  shown  that  in  the  experiments  considered  the  mixing 
efficiency  decreases  rapidly  with  increasing  e/uN2.  Finally,  the  application 
of  our  diffusivity  scaling  to  other  geophysical  flows  is  discussed,  and  it  is 
shown  that  a necessary  requirement  for  the  use  of  this  scaling  is  that  Lt  is 
independent  of  the  rate  of  dissipation  of  turbulent  kinetic  energy. 


Introduction 

A clear  understanding  of  the  irreversible  vertical 
transport  of  mass  in  a stably  stratified  turbulent  flow 
is  fundamental  to  quantifying  the  dynamics  of  density 
stratified  fluids.  The  rate  at  which  this  transport  occurs 
has  been  been  widely  modelled  as  a turbulent  diffusiv- 
ity  for  mass,  Kp.  Since  Kp  influences  the  distribution 
of  heat,  mass,  contaminants  and  biota  throughout  a 
turbulent  fluid  (e.g.  Tennekes  and  Lumley,  1972),  an 
understanding  of  this  quantity  is  essential  to  the  man- 
agement of  aquatic  systems  such  as  lakes,  estuaries  and 
the  oceans. 

The  simplest  model  describing  turbulent  transport 
in  homogeneous  isotropic  turbulence  relates  a turbu- 
lent diffusivity  K to  a turbulent  velocity  scale  U and  a 
turbulent  integral  lengthscale  L ( Tennekes  and  Lumley, 
1972), 

K~UxL.  (1) 


is 


where  u is  the  fluctuating  component  of  velocity  and 
/ is  a linear  dimension  defining  the  scale  of  the  turbu- 
lent field.  If  we  assume  that  U,  u and  L,  l are  well  rep- 
resented by  their  root  mean  square  (ms)  velocity  and 
length  scales  Ut  and  Lt  respectively,  (1)  and  (2)  can  be 
combined  to  give 


(2) 


K ~ e173L*/3.  (3) 

The  rate  at  which  a stratifying  scalar  6 is  irreversibly 
mixed  in  a turbulent  flow  has  generally  been  modelled 
using  an  average  vertical  advective  flux  B'w' , where  9' 
and  w'  are  the  instantaneous  scalar  and  vertical  velocity 
fluctuations  respectively  (see  Gregg , 1987,  for  a full  re- 
view). This  advective  flux  has  been  formulated  in  terms 
of  an  eddy  coefficient  (or  diffusivity)  Kg  as 


6'w' 

Ke  = IT’ 

dz 


This  expression  for  K is  considered  to  represent  a bulk 
eddy  diffusivity.  Taylor  (1935)  found  that  the  rate  of 
dissipation  of  turbulent  kinetic  energy  per  unit  mass,  e, 
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(4) 
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where  88 jdz  is  the  mean  vertical  scalar  gradient  on 
which  the  advective  flux  acts.  Osborn  (1980)  used 
(4)  to  model  Kp  in  the  shear  driven  environment  of 
the  ocean  thermocline.  By  using  the  turbulent  ki- 
netic energy  equation,  Osborn  (1980)  suggested  that 
Kp  < 0.2e/7V2,  where  the  multiplicative  constant  is 
equal  to  Rf/{  1 — Rf),  and  Rf  is  a mixing  efficiency 
taken  to  be  less  than  or  equal  to  0.15.  This  model 
for  Kp  has  often  been  used  in  geophysical  applications 
other  than  originally  intended,  and  with  the  inequality 
written  as  an  equality  such  that 


Kp 


(5) 


Winters  and  D’Asaro  (1996)  argued  that  rather  than 
modelling  Kp  with  (4),  it  can  be  computed  directly  as 


Ka 


<t>d 


Ae. 

dz0 


(6) 


where  (j>d  is  the  irreversible  diffusive  flux  of  density 
across  an  isoscalar  surface  in  a turbulent  flow  and 
dp/dzm  is  the  resorted  density  gradient.  Using  this  defi- 
nition, Barry  et  al.  (2001)  conducted  controlled  labora- 
tory experiments  that  directly  and  independently  mea- 
sured Kp,  the  rate  of  dissipation  of  turbulent  kinetic 
energy  e,  a turbulent  lengthscale  Lt  and  the  buoyancy 
frequency  N2  = —(g/p0)(dp/dz)  in  a closed,  shear  free, 
linearly  salt  stratified  fluid.  Turbulence  was  generated 
in  the  entire  fluid  volume  by  the  steady  motion  of  a 
horizontally  oscillating  vertical  rigid  grid  comprising  1- 
cm  bars  at  a 5-cm  spacing.  A schematic  diagram  of 
the  experimental  configuration  is  shown  in  Figure  1, 
and  full  details  are  given  in  Barry  et  al.  (2001).  Barry 
et  al.  (2001)  described  the  behaviour  of  Kp  with  ref- 
erence to  the  turbulence  intensity  parameter  e/vN2, 
where  v is  the  kinematic  viscosity  of  the  fluid.  Two 
regimes  were  identified:  regime  E (energetic  turbu- 
lence, e/i/N 2 > 300)  and  regime  W (weak  turbulence, 
e/i/N2  < 300).  It  was  found  that 

/f,  = 2to!'V'3(~),/Sand  (7) 

K,=  O.ShW'3^)  (8) 

for  regimes  E and  W respectively.  For  regime  E,  Barry 
et  al.  (2001)  also  found  that 


Figure  1.  Schematic  diagram  of  the  front  view  of  the  ex- 
perimental setup.  The  tank  is  520  mm  long,  and  was  filled 
with  a linearly  salt  stratified  fluid.  The  vertical  rigid  grid 
stirred  the  entire  fluid  volume. 


rms  quantities)  such  that 

Ut  = 3.3  (-jL)  ^ (uN)1'2  Pr~^\  (10) 

where  Pr  = vjn  is  the  Prandtl  number  of  the  fluid. 
This  rms  velocity  scaling  is  the  same  as  that  obtained 
by  using  the  direct  measurements  of  Kp  and  Lt  and 
solving  for  Ut  in  (1).  No  velocity  or  lengthscale  relations 
were  identified  in  regime  W. 

The  primary  focus  of  this  paper  is  to  compare  the 
Barry  et  al.  (2001)  direct  diffusivity  measurements  in 
regime  E with  those  inferred  (via  (4))  using  the  numer- 
ical data  of  Ivey  et  al.  (1998),  the  laboratory  data  of 
Stillinger  et  al.  (1983b),  Rohr  (1985),  and  Itsweire  et 
al.  (1987),  and  the  field  data  of  Saggio  and  Imberger 
(2001).  We  also  use  the  same  data  sets  to  assess  the 
generality  of  the  Barry  et  al.  (2001)  lengthscale  and  ve- 
locity relations  in  the  same  regime. 

The  data  sets 


Lt  = 20 


(vn)1?4 
NU 2 ’ 


(9) 


where  (ok)1/4 /IV1/2  is  the  convective  lengthscale,  Lco. 
This  lengthscale  was  related  to  Ut  using  (2)  (and  assum- 
ing that  l and  u are  well  represented  by  their  respective 


Direct  numerical  simulations 

Using  the  code  described  by  Holt  et  al.  (1992),  Ivey  et 
al.  (1998)  undertook  a series  of  direct  numerical  simula- 
tions (DNS)  of  stratified  shear  flows  using  a 1283  point 
computational  domain  in  which  N2  and  the  square  of 
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the  mean  vertical  shear  S2  = ( dU/dz )2  were  held  con- 
stant. Although  only  one  value  of  S2  was  used  for 
all  the  simulations  discussed  here,  N2  was  varied  be- 
tween simulations  such  that  the  Richardson  number, 
Ri  =s  N2/S2,  varied  from  0.075  to  1.0.  The  experiments 
were  run  for  non-dimensional  shear  time,  0 < St  < 8, 
however  we  consider  the  results  from  only  St  > 2 where 
the  turbulence  becomes  self-adjusted  and  has  lost  mem- 
ory of  the  initialisation  ( Ivey  et  al.,  1998).  At  each  time- 
step  within  a simulation,  all  turbulent  quantities  were 
calculated  as  ensemble  averages  over  the  entire  fluid  vol- 
ume. With  the  application  to  geophysical  flows  in  mind, 
we  consider  only  the  DNS  experiments  with  Pr  > 1. 

Laboratory  experiments 

The  experiments  of  Stillinger  et  al.  (1983b),  Rohr 
(1985),  and  Itsweire  et  al.  (1987)  measured  turbulent 
quantities  at  different  points  downstream  of  a fixed  grid 
at  the  head  of  a water  tunnel  (see  Stillinger  et  al .,  1983a, 
for  a full  description  of  the  experimental  configuration). 
These  experiments  used  a variety  of  grids  with  bars 
of  diameter  0.318  to  0.635  cm  at  spacings  of  1.905  to 
3.81  cm.  In  all  the  experiments,  salt  stratified  water 
was  continually  recirculated  through  the  tunnel,  allow- 
ing the  collection  of  timeseries  data  at  several  positions 
downstream  of  the  stationary  grid.  Average  turbulent 
quantities  were  computed  at  these  positions  by  aver- 
aging the  collected  timeseries  at  each  point.  The  Rohr 
(1985)  experiments  included  mean  shear.  From  numeri- 
cal simulations,  Itsweire  et  al.  (1993)  found  that  in  these 
experiments,  e was  underestimated  by  approximately  a 
factor  of  2.  We  have  applied  this  correction  factor  to 
the  data  sets  for  this  analysis. 

Field  experiments 

Saggio  and  Imherger  (2001)  measured  turbulent  prop- 
erties in  the  metalimnion  of  Lake  Kinneret,  Israel,  dur- 
ing three  consecutive  summer  periods.  During  these 
experiments,  a portable  flux  profiler  (PFP)  capable  of 
resolving  all  three  turbulent  velocity  components  was 
traversed  vertically  through  the  water  column  at  ap- 
proximately 10  cm/s.  This  probe  also  allowed  the  mea- 
surement of  the  instantaneous  vertical  advective  buoy- 
ancy flux  p'w' , and  the  rate  of  dissipation  of  turbulent 
kinetic  energy,  e. 

Results 

Turbulent  lengthscales  and  velocities 

In  Figure  2 we  plot  Ut  against  (eLt)1/3,  for  the  above 
numerical  and  laboratory  data  sets  where  e/vN2  > 300. 
This  rms  quantity  was  computed  as 


Figure  2.  Root  mean  square  turbulent  velocity  Ut  vs 
(eLt)1/3,  in  cm/s;  • , Stillinger  et  al.  (1983b);  o , Rohr 
(1985);  + , Itsweire  et  al.  (1987);  * , Ivey  et  al.  (1998). 
The  solid  line  is  slope  1,  with  the  equation  given  by  (13). 


Ut  = V2(u'2)  + w'2  (11) 

for  the  laboratory  data  and 

Ut  = vV2  + v'2  + w'2  (12) 


for  the  DNS  data.  The  solid  line  in  Figure  2 is  slope  1 
and  has 

Ut  = 1.5  (eLt)1/3  (13) 

consistent  with  Ivey  et  al.  (1998).  Figure  3 shows  the 
normalised  turbulent  velocity  scale  against  e/vN2  for 
the  same  data,  where  the  solid  line  is  given  by 


Ut 

(vN)U2Pr~1/12 


(14) 


Equations  (13)  and  (14)  imply  that  the  rms  turbulent 
lengthscale  Lt  in  the  DNS  and  laboratory  data  is  de- 
scribed by 


Lt  = 35 


(i/k)1/4 

TV1/2 


— 35  Lco. 


(15) 


Similarly,  Saggio  and  Imberger  (2001)  found  for  field 
data  in  a heat  stratified  lake  that 


Ut 

{vN)l/2Pr~l/ 12 


(16) 


and 

Lt  = 15  Lco,  (17) 

where  Ut  was  computed  using  (12). 
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Figure  3.  Normalised  turbulent  velocity 

Ut/((vN)1/2Pr~l/12)  vs  e/uN2;  ■ , Stillinger  et  al. 

(1983b);  o , Rohr  (1985);  + , Itsweire  et  al.  (1987);  * , Ivey 
et  al.  (1998).  The  solid  line  is  slope  1/3,  with  the  equation 
given  by  (14). 


Figure  4.  {b/N2)/(v2^K1^)  vs  e/vN2;  ■ , Stillinger  et  al. 
(1983b);  o , Rohr  (1985);  + , Itsweire  et  al.  (1987);  * , Ivey 
et  al.  (1998).  The  equations  of  the  lines  are  given  in  the 
text. 


Buoyancy  flux  estimates  of  Kp 

In  Figure  4 we  plot  the  quantity  ( b/N2)/{v2^3n 1^3) 
against  e/vN2,  where  b = (g/ p0)p'w'  is  the  buoyancy 
flux.  The  solid  lines  are  given  by 


....  b — 24  (— Y/3 

7V2(i/2/3«1/3)  \vN2) 

for  the  laboratory  data  and 

b y/3 

N2^2/3k1/3)  °Vl/lV2/ 


(18) 


(19) 


for  the  DNS  data  where  e/i jN2  > 300. 

Saggio  and  Imberger  (2001)  found  from  field  data  in 
Lake  Kinneret  that  when  e/uN2  > 36, 


N2^2/3k1/3) 


/ e \ V3 

- 12  (^)  ■ <2°) 


Discussion 

The  functional  forms  of  equations  (14),  (16)  and  (15), 
(17)  are  consistent  with  the  Barry  et  al.  (2001)  scalings 
in  (10)  and  (9),  respectively.  In  particular,  (15)  and  (17) 
demonstrate  that  in  all  these  data  sets  the  turbulent 
lengthscale  is  not  a function  of  the  rate  of  dissipation  of 


turbulent  kinetic  energy  when  e/vN2  > 300.  Further, 
these  velocity  scales  are  consistent  with  the  results  of 
Barry  et  al.  (2001)  who  suggested  that  for  a given  N in 
regime  E,  an  increase  in  e does  not  result  in  an  increase 
in  the  turbulent  overturn  scale  but,  rather,  leads  to  an 
increase  in  the  magnitude  of  the  rms  turbulent  velocity. 

Although  we  note  some  variation  in  the  constant 
coefficient,  the  agreement  between  the  turbulent  ve- 
locity and  lengthscale  relations  of  Barry  et  al.  (2001) 
and  those  of  the  data  sets  presented  here  suggests  that, 
to  first  order,  the  configuration  of  the  turbulence  gen- 
eration mechanism  does  not  influence  the  physics  of 
the  turbulence  when  e/uN2  > 300.  Since  all  data 
sets  are  reasonably  well  described  by  Lt  ~ Lco  and 
Ut/((vN)'/2Pr- J/12)  ~ (e/rdV2)1/3,  we  assume  that 
(1)  also  applies  to  the  experiments  detailed  above  and, 
if  it  were  directly  measured,  Kp  would  be  well  described 
by  (7)  in  regime  E. 

Figure  4 shows  that  for  regime  E there  is  good 
agreement  between  (7)  (given  by  the  solid  line)  and 
b/(N2i/2^3K1^3)  computed  using  the  laboratory  data. 
Although  there  is  some  scatter  in  the  laboratory  data 
(consistent  with  the  scatter  in  the  velocity  and  length- 
scale  measurements  shown  in  Figures  2 and  3),  the  trend 
in  b/{N2v2/3K 1/3)  is  consistent  with  (7)  in  regime  E.  We 
conclude  that  in  these  laboratory  experiments,  (4)  can 
be  used  to  satisfactorily  approximate  the  true  diffusiv- 
ity  measured  by  Barry  et  al.  (2001)  when  e/vN2  > 300. 
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Figure  5.  (6/iV2)/(i/2/,3/c1^3)  vs  e/uN2;  ■ , Stillinger  et  al. 
(1983b);  o , Rohr  (1985);  4-  , Itsweire  et  al.  (1987).  The 
solid  line  is  (0.2e/N2)/(i/2/V/3),  with  Pr  = 700. 


The  quantity  b/(N2v2l3K1/3)  computed  from  both 
DNS  and  Saggio  and  Imberger  (2001)  data  shows  the 
same  functional  dependence  on  (e/uN2)1/3  as  (7)  when 
e/uN2  > 300.  This  implies  that  the  physics  of  the  tur- 
bulence in  regime  E (i.e.  the  relationship  between  Kp , 
Ut  and  Lt ) is  well  described  by  the  scale  forms  of  the 
Barry  et  al.  (2001)  results.  We  note,  however,  that  the 
constant  coefficients  in  the  DNS  and  field  data  sets  are 
different  from  those  of  the  Barry  et  al.  (2001)  measure- 
ments. It  is  currently  not  clear  why  this  should  be  the 
case. 

Figure  5 shows  the  estimates  of  Kp  from  the  labora- 
tory data  sets  using  (4)  compared  with  the  prediction 
from  Osborn  (1980)  in  (5)  for  Pr  = 700.  The  figure 
shows  that  for  e/uN 2 > 300,  the  model  predictions 
from  (4)  and  (5)  diverge.  For  highly  energetic  turbu- 
lence, the  two  estimates  are  widely  disparate,  such  that 
at  ejvN 2 « 105  they  are  different  by  more  than  two 
orders  of  magnitude. 

Barry  et  al.  (2001)  argued  that  the  most  likely  cause 
of  this  divergence  for  large  ejvN2,  is  the  selection  of  a 
constant  mixing  efficiency  Rf  in  deriving  (5).  Ivey  and 
Imberger  (1991)  described  the  behaviour  of  this  mixing 
efficiency  in  terms  of  the  turbulent  Froude  number, 

■ (21> 


and  the  turbulent  Reynolds  number, 

<22> 

where  L0  — (e/N3)1^2  is  the  Ozmidov  scale,  Lc  is  a 
centred  displacement  scale  equivalent  to  Lt,  and  Lk  = 
(^3/e)1/4  is  the  Kolmogorov  scale.  In  particular,  Ivey 
and  Imberger  (1991)  found  that  Rf  has  a maximum 
of  approximately  0.2  at  Frx  ~ 1.2  (or  La  « Lt),  and 
that  Rf  decreases  rapidly  when  Frx  deviates  from  this 
value,  such  that  for  Frx  > 1.2, 

R>  - irk*  (23) 

for  large  Rex-  If  we  assume  that  Lt  (and  hence  Lc)  is 
well  described  by  27.5 Lco  for  e/uN2  > 300  (where  the 
coefficient  of  27.5  is  a representative  value  from  the  Still- 
inger et  al.  (1983b),  Rohr  (1985),  Itsweire  et  al.  (1987) 
and  Barry  et  al.  (2001)  data  sets),  we  can  substitute 
27. 5 L co  for  Lc  in  (21)  and  (22)  to  give 

FrT  = 0.11  Prl/ 2)  and  (24) 

ReT  = S3{^Pr  ) • <25) 

respectively.  Ivey  and  Imberger  (1991)  argued  that  any 
two  of  the  non-dimensional  parameters  FrT,  Rex  and 
e/uN2  are  necessary  and  sufficient  to  describe  the  dy- 
namics of  a general  turbulent  field.  However,  equations 
(24)  and  (25)  show  that  for  regime  E in  the  above  ex- 
periments we  need  only  one. 

We  plot  Frx  against  e/vN2Pr1^2  in  Figure  6,  for 
e/t/N2  > 300,  where  the  solid  line  is  given  by  (24). 
Equations  (23)  and  (24)  show  that  Rf  necessarily  de- 
creases with  increasing  ejvN2, 


Rf  — 


1 

l + 0.04(^)2/3Pr1/3' 


(26) 


Figure  7 shows  Rf  against  e/vN2  for  these  data  sets, 
and  clearly  demonstrates  that  a diffusivity  modelled  as 
Kp  = Rf  /(1—Rf)e/N2  with  Rf  constant  at  0.15  (shown 
in  Figure  7 as  a horizontal  solid  line)  will  overestimate 
the  true  diffusivity  at  large  ejvN2  in  these  experiments. 
This  result  is  consistent  with  the  divergence  of  the  es- 
timates of  Kp  from  (4)  and  (5)  shown  in  Figure  5. 


Application  to  geophysical  flows 

Equations  (24)  and  (26)  hold  if,  and  only  if,  the  tur- 
bulent lengthscale  Lx  is  well  described  by  Lco.  Sev- 
eral field  investigations  however,  (e.g.  Crawford,  1986; 
Gregg , 1987;  Peters  et  al,  1988;  Imberger  and  Ivey, 
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Figure  6.  Turbulent  Froude  number  Frr  vs 
( e/vN2Pr 1/2)1/3;  • , Stillinger  et  al.  (1983b);  o , Rohr 
(1985);  + , Itsweire  et  al.  (1987);  A , Barry  et  al.  (2001);  x , 
Ferron  et  al.  (1998).  The  solid  line  is  0.11(«/WV2Pr1/2)1/3. 


Figure  7.  Rf  vs  e/vN2;  ■ , Stillinger  et  al.  (1983b);  o , Rohr 
(1985);  + , Itsweire  et  al.  (1987);  A , Barry  et  al.  (2001). 


1991;  Peters  et  al.,  1995;  Ferron  et  al.,  1998)  have  found 
that  in  some  geophysical  flows  Lt  ~ L0,  or  equivalently, 
that  Frr  « 1.  In  these  situations,  this  relationship  be- 
tween Lt  and  L0  has  been  found  to  hold  regardless  of 
the  value  of  e/vN2.  For  instance,  Ferron  et  al.  (1998) 
observed  that  in  the  Romanche  fracture  zone, 

L„  = 0.95(±0.6)£t,  (27) 

even  when  e/vN 2 ss  500,000,  where  Lt  is  the  Thorpe 
displacement  scale  and  is  a measure  of  Lt-  The  corre- 
sponding value  of  Frx  is  order  1 and  is  shown  in  Fig- 
ure 6 as  an  x.  The  Ferron  et  al.  (1998)  measurements 
clearly  do  not  conform  to  the  trends  of  the  data  sets 
previously  considered  here,  fundamentally  because  Lt 
does  not  scale  like  Lco. 

The  measurements  of  Ferron  et  al.  (1998)  imply  that 
Frr  is  always  approximately  unity,  regardless  of  the 
magnitude  of  e/vN2.  The  arguments  of  Ivey  and  Im- 
berger  (1991)  then  require  that  the  mixing  efficiency  in 
the  Romanche  fracture  zone  is  constant  and  maximal, 
and  hence  the  large  value  of  e/vN2  cannot  be  used  in 
isolation  to  infer  that  Rf  is  vanishingly  small.  Rather, 
either  FrT  or  Rer  must  also  be  included  in  a description 
of  the  turbulent  dynamics.  In  such  cases,  the  scaling  re- 
sults of  Barry  et  al.  (2001)  should  not  be  expected  to 
hold. 

Conclusions 

We  have  compared  the  turbulent  lengthscale,  ve- 
locity and  diffusivity  scalings  presented  by  Barry  et 
al.  (2001)  with  data  from  Ivey  et  al.  (1998)  DNS  ex- 
periments, laboratory  experiments  of  Stillinger  et  al. 
(1983b),  Rohr  (1985),  and  Itsweire  et  al.  (1987)  and 
the  field  experiments  of  Saggio  and  Imberger  (2001). 
For  e/vN 2 > 300,  the  Barry  et  al.  (2001)  velocity  and 
lengthscale  results  are  good  descriptors  of  the  physics 
of  the  turbulence  in  all  these  experiments,  to  constant 
coefficient.  We  have  shown  that  the  turbulent  length- 
scales  in  all  these  experiments  in  this  regime  are  inde- 
pendent of  the  rate  of  dissipation  of  turbulent  kinetic 
energy.  This  in  turn  suggests  that  for  a given  N,  an  in- 
crease in  e does  not  increase  the  turbulent  lengthscale, 
but  rather,  leads  to  an  increase  in  the  magnitude  of  the 
turbulent  velocity. 

We  have  shown  that  when  e/vN2  > 300,  a diffusivity 
modelled  using  (4),  is  a good  representation  of  the  direct 
measurements  of  Barry  et  al.  (2001)  for  the  experiments 
considered,  to  constant  coefficient. 

We  have  discussed  the  parameterisation  of  (4)  as 
Kp  — 0.2 e/N2  and  have  shown  that,  in  the  data  sets 
considered,  this  overestimates  the  true  diffusivity  for 
e/vN2  > 300.  We  have  suggested  that  this  is  a conse- 
quence of  selecting  a constant  mixing  efficiency  Rf  = 
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0.15,  as  for  increasing  e/vN 2,  Rf  necessarily  decreases 
in  the  data  sets  considered. 

Finally,  we  have  discussed  the  application  of  the 
Barry  et  al.  (2001)  scalings  to  geophysical  flows  and 
have  suggested  that  a necessary  requirement  for  the 
use  of  these  scalings  is  that  the  turbulent  lengthscale 
Lt  is  independent  of  e.  If  this  constraint  is  satisfied 
(i.e.  Lt  ~ Lco)  then  Kp,  and  indeed  all  the  turbulent 
dynamics,  can  be  described  in  terms  of  only  a single  pa- 
rameter e/vN2 . When  this  lengthscale  constraint  is  not 
satisfied,  the  Barry  et  al.  (2001)  scalings  do  not  apply. 
Future  work  must  now  investigate  the  circumstances 
under  which  this  lengthscale  constraint  is  satisfied  in  a 
variety  of  geophysical  flows. 
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